Comment on "Reverse Brazil Nut Problem: Competition between 
Percolation and Condensation" 

In a recent paper j[J the Brazil nut problem (BNP), where particles with large diameters 
rise to the top when exposed to vertical shaking and the reverse Brasil nut problem (RBNP) 
where the large particles segregate to the bottom were investigated. In realistic situations 
these effects are driven or at least accompanied by global convection. Convection rolls are 
observed already in mono-disperse systems: appearantly vibrations lead to temperature 
gradients with temperature defined through the average kinetic energy of the particles. In 
order to suppress convection the authors present Molecular Dynamics (MD) simulations for 
binary hard sphere mixtures at fixed temperature. It appears legitimate to first consider 
this simplified situation although it might be rather remote from reality. Depending on the 
mass and diameter ratios of the two hard sphere species, the simulations show pronounced 
BNP and RBNP behaviours. For an explanation a new condensation versus percolation 
driven segregation mechanism was put forward, which has attracted much attention in the 
press [0. The purpose of this comment is to show that the proposed mechanism is incorrect. 
All the effects observed in the MD simulations can be understood by simple conventional 
thermodynamics. 

Because in the solid phase particles may no longer exchange positions |J, we conjecture 
that demixing must happen already in the fluid phase. Therefore we consider a binary hard 
sphere mixture with species A and B having masses vha > tub and diameters > ds in 
the liquid phase with an appropriate equation of state (EOS), e.g. the Mansoori-Carnahan- 
Starling-Leland equation |J, with £j = (n^/n) d\ + (ns/n) d B 
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where n = ua + n B is the particle density and x — xa + xb = {ir/&)d\riA + (ir/6)dnnB the 
packing fraction. It is then straightforward to compute the Helmholtz free energy Q under 
gravity 

F = E-T-S. (2) 

Variation of F with respect to the particle densities ua and ns for fixed particle numbers 
Na and Nb leads to two coupled nonlinear equations which are solved numerically. In the 
limit NAd\ <§C Nsd B and dA 3> ds these equations respect Archimedes' principle. Note 
however, that the initial condition of equal layer numbers N^d^ = Nsdg used in Q] is 
incompatible with this limit. Consequently Archimedes' principle need not be obeyed in the 
presented MD simulations in contrast to our naive expectation. Calculated filling fractions 
Xa{z) and xb{z) as functions of the height z are depicted in Fig. 1 for various mass ratios 
mA/mB keeping the diameter ratio d^/ds — 2 and particle number ratio Na/Nb = 1/4 
fixed. The density profiles compare very well with the MD simulations presented in |Q. 
For mA/rriB = 2 we encounter the BNP with almost pure B-liquid at the bottom. The 
mA/nriB = 6 case corresponds to the RBNP with an incomplete demixing at the bottom. 
For mA/nris = 4 the two liquids remain mixed. All these details agree with the corresponding 
MD simulations. But note that in our calculation the temperatures T > Tc > Ta lie above 
the temperature Tc where the binary liquid begins to solidify: the total packing fraction 
lies sufficiently below close packing everywhere. In contrast the temperatures used in Q are 
supposed to be quenched between Tb < T < Ta {Ta and Tb represent the temperatures 
where condensation sets in for the separated species || , see below) such that parts of the 
sample are " frozen" . 
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Figure 1: Calculated filling fractions xa and xb of binary hard sphere liquids under gravity 
for various mass ratios. The length unit equals the container size used in the MD simulations 
of Ref. Q. The BNP and RBNP are clearly recognized. With decreasing temperature the 
profiles are shifted to smaller hights z. 
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We must conclude that demixing indeed occurs in the liquid phase. It is clearly due 
to the competition between gravity and entropy rather than to that of condensation and 
percolation. There is nothing mysterious or special going on and all the effects are understood 
by conventional thermodynamics. Finally, if the temperature is further decreased, the binary 
liquid eventually begins to condensate locally according to its filling fractions xa and 15 in 
the fluid phase. In the RBNP this starts from the bottom whereas in the BNP, depending 
on the parameters, species A in the middle may condensate first (cf. Fig. 1). However, 
there exists no condensation driven demixing. 

In order to discuss the transition from the BNP to the RBNP and also to understand the 
effects on an even simpler basis, we may consider the two species separated: A above B for 
the BNP and B above A for the RBNP. The advantage is that we may then use the much 
simpler EOS for one-component hard sphere fluids e.g. the Carnahan-Starling expression 



P = nTZ(x), Z(x)= 1 + X {i + _ x2 x)3 X " (3) 

obtained from (I) by setting dA — ds- The details of the EOS are unimportant here, in 
fact, every Z(x) which grows sufficiently fast with x (eigen-volume effect) leads to similar 
results. We assume that at a filling fraction xo ~ 0.637 equal to random close packing [Q, 
instead of FCC or HCP close packing xq = n/>/l8 — 0.74048, condensation sets in, which 
leads to temperatures 
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(gravity g) which exceed those in || by roughly a factor of 3. We choose xa(0) — xo at the 
bottom of the RBNP which corresponds to the temperature Trbnp = Ta + (dA/ds) 3 Tb > 
Ta+Tb where this system is just still fluid. Then also the BNP configuration is fluid, because 
there " freezing" sets in at lower temperatures (either at Tb np = Ta at the boundary of the 
two species or at Tbnp = (ds/cU) 3 Ta + Tb at the bottom, depending on the parameters 
used). With the above choice the filling fractions at the bottom and at the boundaries are 
fixed by the corresponding pressures and the Helmholtz free energies Fbnp and Frbnp 
may be computed. For Fbnp < Frbnp the free energy favors the BNP and vice versa. 
We consider the curves Fbnp = Frbnp which separate BNP from RBNP for the following 
cases: 

(i) equal number of particles Na = Nb 

(ii) equal number of layers NaA\ — Nsd B 

(iii) equal volumes NAd\ = Nsd B 

(iv) equal masses Nattla — Nsms- 

The last case is not explicitly shown in Fig. 2 because it almost coincides with case (ii). The 
MD simulation results, case (ii), are also shown with circles (BNP), boxes (RBNP) and 
triangles (mixed). In the vicinity of the curves mixing is expected, demixing occcurs for 
parameter combinations lying sufficiently far away. We find that the particle numbers are 
of minor influence. This is expected: the BNP and RBNP should not be easily reversed by 
adding or subtracting a few particles of either species. 
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Figure 2: Curves which separate BNP from RBNP for different particle number ratios in 
the mass ratio and diameter ratio parameter space. The MD simulation results of Rcf. |IJ 
are also shown. 

Lets contrast our results with the "percolation-condensation" formula 

Ta = N A m A d\ = d% 

T B N B m B d% d% 1 ' 

of Ref. (!]]. This formula seems to give reasonable results dA/d B = (mi/ms) 1 / 2 for the 
equal-number-of-layers case (ii), but certainly a straight line would not be worse. For equal 
volumes, case (hi), (5) leads to d A /d B — (m/i/m^) 1 / 3 and the mA/rn B ~ 6, d A /d B = 2 
system is predicted BNP, which is incorrect. The calculated profiles look very similar to 
the corresponding ones of case (ii) (Fig. 1) with the only difference that x B is somewhat 
larger because of the relatively larger particle number N B . An MD calculation certainly 
would confirm this. For equal particle numbers, case (i), the BNP is obtained for itia < m B 
wheras the equal mass case (vi) is predicted to be always mixed. All these implications are 
wrong and we must conclude that the formula (5) works only approximately for case (ii) by 
pure accident. 

Preparing this comment, I realized that J. A. Both together with one of the authors of 
Ref. [0, D.C. Hong, have just submitted a paper to the archive ||], which treats the problem 
along similar lines. This paper, however, lacks the clear statement that the "condensation 
- percolation theory" is incorrect. In view of the publicity given to this picture Q such a 
statement would seem appropriate. 

H. Walliser 

Fachbereich Physik, Universttat Siegen, D-57068 Siegen, Germany 
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